Nagel Scaling and Relaxation in the Kinetic Ising Model on a n- Isotopic Chain 
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The kinetic Ising model on a n-isotopic chain is considered in the framework of Glauber dynamics. 
The chain is composed of N segments with n sites, each one occupied by a different isotope. Due 
to the isotopic mass difference, the n spins in each segment have different relaxation times in the 
absence of the interactions, and consequently the dynamics of the system is governed by multiple 
relaxation mechanisms. The solution is obtained in closed form for arbitrary n, by reducing the 
C/3 ' problem to a set of n coupled equations, and it is shown rigorously that the critical exponent z 

is equal to 2. Explicit results are obtained numerically for any temperature and it is also shown 
that the dynamic susceptibility satisfies the new scaling (Nagel scaling) proposed for glass-forming 
liquids. This is in agreement with our recent results (L. L. Gongalves, M. Lopez de Haro, J. Tagiiena- 
I , Martinez and R. B. Stinchcombe, Phys. Rev. Lett. 84, 1507 (2000)), which relate this new scaling 

function to multiple relaxation processes. 
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The search for universal behavior in various physical model systems has been one of the clues to uncover fundamental 
regularities of nature. About ten years ago, a scaling hypothesis (the Nagel scaling) was proposed to describe 
experimental work on dielectric relaxation in glass- forming liquids m. This hypothesis was meant to replace the 
more usual normalized Debye scaling, where the frequency is scaled with the inverse of the (single) relaxation time 
and the real and imaginary parts of the dynamic susceptibility are then divided by their values at zero and one, 
respectively. However, and in spite of its phenomenological success and apparent ubiquity, its physical origin is as yet 
. not quite well understood. 

Very recently, we put forward what we believe to be the first bona fide microscopic model in which the Nagel 
scaling is shown to arise In this previous work, we computed the magnetization, the dynamic critical exponent 
z and the frequency and wavevector dependent susceptibility of the kinetic Ising model on an alternating isotopic 
chain with Glauber dynamics || . Our results indicated that as soon as the two relaxation times of this model became 
substantially different, the agreement with the Nagel scaling improved significantly. The question remained of whether 
the inclusion of more relaxation mechanisms would lead to the same sort of results. It is this issue that we mainly 
want to examine in this paper. Hence we present an extension of the alternating isotopic chain in which rather than 
q ' only two relaxation times, it is n relaxation times which are involved. 

The model consists of a linear chain with N segments, each one containing n sites occupied by isotopes characterized 
by n different spin relaxation times. The Hamiltonian for the I th segment is the usual Ising Hamiltonian given by 



X 

where aij is a stochastic (time-dependent) spin variable assuming the values ±1 and J the coupling constant. Note 
that if n = 2, this model reduces to the alternating isotopic chain treated in Ref. The configuration of the 
segment is specified by the set of values {07,1, 07,2, ...<Ji ;n \ = |<7 !,n } at time t. This configuration evolves in time due 
to interactions with a heat bath. We assume for the segments the usual Glauber dynamics so that the transition 
probabilities are given by 

= \ai (l - l^i-ivi, + Wi+x)) , (2) 
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where 7 = tanh (2 J/fc^T), ks being the Boltzmann constant and T the absolute temperature, and is the inverse 
of the relaxation time Tj of spin i in the absence of spin interactions. 

The time dependent probability P ({ct''™} , t) for a given spin configuration satisfies the master equation 



dP({a l ' n ) ,t) 



dt 
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+ J2 w u (-a M )P(T H {^"},i), 



(3) 



where TJ, jcr''™} = 07.2, 07,1+1, ■■■0'i,n}- The dynamical properties we are interested in, namely the 

magnetization, the dynamic critical exponent and the susceptibility, require the knowledge of some moments of the 
probability P ({ct''™} , t). Hence, we introduce the following expectation values defined as: 



qu (t) = {a u (t)> = J2 <n,iP {W' n } ,t) 



(4) 



where the sum runs over all possible configurations. Using the master equation (3) and this definition of qu (t), one 
can easily derive the result 



dqi. 



dt 
dqi,j 

dt 
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dt 
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)) , (2 < j < n - 1) 



Introducing now the Fourier transform qQj (I — 1, 2, ...n) defined by 



1 N 

q Q ,i = —= V enp(iQld)q, 



where Q = (m = 0, 1, 2...N), d — na 7 a being a lattice parameter and the vector given by 



one can derive the following result 
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where the matrix Mq has a rather suggestive structure given by 
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The solution to Eq. (8), which yields the magnetization, is straightforward, namely 

* Q (t) = e M «*v!/ Q (0) . 



(5) 



(6) 
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(10) 
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The relaxation process of the wave- vector dependent magnetization is determined by the eigenvalues of Mq denoted 
as Xiq (I — 1, 2, ...n) and which for n > 5 have to be computed numerically. The same applies to the dynamic critical 
exponent z, which is obtained by imposing the scaling relation tq ~ £, z f(£,Q) for the critical mode, Xiq, with the 
(Q-dependent ) relaxation time tq — ™xj— > m ^ ne re gi° n T ^ and Q — > where Xiq — > 0, and the correlation 
length £ is £ ~ exp(2J/k B T). This yields 

where C is an irrelevant constant. By plotting In ^— -^-^j vs. and considering the limit T — > 0, we have checked 

analytically for n = 2 and numerically for n — 3,4 and 5 and various values of the a's that z = 2, so that this 
model belongs to the same universality class of the uniform Ising chain. 

Let us now introduce the spatial Fourier transform cq(£', t' + t) of the time-dependent correlation defined by 



N N n n 

c Q (t', *' + *) = ^£EEE e- iQdl e^' (a^ (f ) a V j (f + t)) . (12) 

1=1 V = l i=l j=l 

Then, the t' — > oo limit of the temporal Fourier transform of cq it' , t' + t) , denoted by Cq(u>), is given by 
^ ^ 7i n „oo 

C Q H = Hm — E E / + *)> exp(-iojt)dt, (13) 

with OQ^m — -j= ^2 ex p(*Q^) cr ;,m- After some lengthy and tedious but not very difficult algebra, using eq. (Q) and 
v i=i 

the above definitions one may derive the result 
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where 



n n n 



91Q = E E E (°-Q,i°Q,i)e q ■ ( 15 ) 
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Here, the static correlation functions (c-q^ctq^) e are given by 
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while u = tanh(^j^r), a m ; denotes the (to, — element of the matrix A formed with the eigenvectors of Mq and ay 
the (l,j)— element of the matrix A" 1 . Finally, by using the fluctuation dissipation theorem the response function 
Sq (uj) turns out to be given by 



1 / 1 - - 
Sq (w) = jrj; -EE^ft 1 ^ - iwC Q^) 
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It should be noted that the frequency dependent susceptibility is x ( w ) = kBTSo(ui)(l — u)/(l + u). Thus, 

x(w) = Hziy; T ^_ (is) 

u + 1 ^— ' ioj — Xiq K ' 

i=i u 



3 



If we set all the a's equal in Eq. (17), i.e. we take the uniform chain, then of course the resulting susceptibility 
has the simple Debye form. In any case, its general structure is a linear combination of n Debye-like terms. In order 
to investigate whether the Nagel scaling holds for this susceptibility, it is convenient to recall that in the so-called 
Nagel plot the abscissa is (1 + W)log 10 (u/ui p ) /W 2 and the ordinate is log 10 (x"(u) w p/ w/ ^x) /W. Here, x" is the 
imaginary part of the susceptibility, W is the full width at half maximum of x", uj p is the frequency corresponding 
to the peak in x", and A\ = x(0) — Xoo is the static susceptibility. For the sake of illustration in Figs. 1 and 2 we 
present Nagel plots for the cases n = 2 (with ai = 1 and 012 = 2), and n — 3 (with oti = 1 and 0:2 = 2 and az = 3), 
respectively, and different values of the reduced (dimensionless) temperature T* = ksT /2J . By comparing the two 
figures, we can see that even for this case, where the relaxation times are very close, there is a marked improvement 
of the scaling in the low temperature limit with the increase of the relaxation mechanisms. On the other hand, the 
Debye scaling, shown as inset in the figures, presents an opposite behaviour. This result gives further support to the 
hypothesis that the Nagel scaling is related to multiple relaxation mechanisms, as discussed in our previous work [2]. 
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Figure captions 

Figure 1. Nagel plot for n = 2 (with ct\ = 1 and ct2 = 2) and for T* — 5, 7, 10, 50, 100 and 150. There is reasonable 
agreement with the scaling form for this choice except for low T* . The insert contains the plot of x" (^) / x" (^p) vs - 
ui/uip in order to test the Debye-like behavior. 

Figure 2. The same as Fig. 1 but with the choice n = 3 (ax = 1 , a 2 = 2 and a 3 = 3) and T* = 5, 7, 10, 50, 100 
and 150. The improvement in the agreement with the Nagel scaling is rather noticeable, while the opposite trend is 
observed with respect to the Debye scaling. 
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